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Abstract 

Convective contamination of self-diffusion experiments with an applied magnetic field is considered 
using a two-dimensional axisymmetric model. Constant, uniform, and an additional non-unifonn heat 
fluxes are imposed along the sidewall of the cylinder while constant heat loss occurs through the top and 
bottom. In this model, due to a very small thermal Peclet number, convective heat transfer is neglected, 
and the flow is steady and inertialess. Time-dependent concentration is solved for various values of the 
mass Peclet number, Pe , (the ratio between the convective transport rate and the diffusive transport rate) 
and different magnetic field strengths represented by the Hartmann number Ha. Normalized values of 
these diffusivities vs. effective Pe are presented for different imposed temperature profiles. In all cases, 
the diffusivity value obtained through the simulated measurement increases as the effective Pe m increases. 
The numerical results suggest that an additional periodic flux, or “hot” and “cold” spots, can significantly 
decrease the convective contamination in our geometry. 

1. Introduction 

During liquid self-diffusion experiments, a small non-uniformity in temperature within the melt may 
drive buoyant convection. This convection can result in erroneous values of the measured diffusivity. 
Verhoeven [1] emphasized that any horizontal component of a density gradient in the liquid results in 
spontaneous convection with no threshold. Alexander et al. [2] showed that for three-dimensional (3D) 
time-dependent transport in the presence of gravity with horizontal temperature non-uniformities across 
the sample as low as 1 and 0.1 K, convective transport rates in 1 and 3 mm capillaries, respectively, can 
exceed diffusive transport rates, and, thus, result in higher values of measured diffusivities. Alexander and 
Banish [3] presented the results of combined numerical modeling and order-of-magnitude estimates of the 
sensitivity of convective contamination to low-gravity in self-diffusion experiments. 

Applied magnetic fields have been used to suppress the buoyant convection in the liquid metals and 
semiconductors due to their large electrical conductivities [4-6]. Youdelis et al. [7] showed that at high 
magnetic field strengths, the diffusion process itself can be significantly modified through the Lorentz 
body force acting on the conducting liquids ions and electrons. Alboussiere et al. [8] found that with 
electromagnetic damping the convective contribution to the effective mass Peclet number is scaled as 
Hartmann number Ha' 4 . 
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In actual self-diffusivity experiments, measured diffusivity results that deviate by less than 5% from the 
(known) real value, D (p would be considered “acceptable. The purpose of this paper is to quantitatively 
estimate the magnitude of allowable temperature non-uniformities in the liquid that will guarantee that 
measured diffusivity values lie within 5% of the actual value when the experimental system is subject 
to an axial magnetic field. To model the experiment, axisymmetric, time-dependent transport due to a 
combination of steady flow and diffusion is considered. In this basic model, a steady, uniform heat flux 
that produces a radial temperature difference in the liquid is imposed on the sidewall of the cylindrical 
container while unifonn heat losses are imposed on the top and bottom walls. Then, a steady, spatially- 
periodic heat flux is superimposed on the sidewall of the cylinder to predict the effect of spatial temperature 
variations which represent the localized hot and cold spots along the sidewall. Note that the overall flux 
for this additional periodic condition is zero. Simulated diffusivity results for two different magnetic field 
strengths and five various temperature profiles along the sidewall are presented in this paper. 

2, Problem Formulation 

In this model, a Boussinesq fluid contained in a closed vertical circular cylinder of length Z = 30 mm 
with an inside radius R of 1.5 mm. Gravity acts downward along the cylinder axis while a unifonn axial 
magnetic field is applied in the opposite direction as in Fig. la. The origin lies at ( R = 0) and z = Z/(2R). 
The isopicnic radioactive tracer is initially located at the bottom of the cylinder. The dimensionless model 
has the top and bottom limits of z =10 and -10, respectively, while the vertical wall lies at r = 1 . 


2.1. Thermal Problem 

For sufficiently large values of magnetic field strength B , the magnetic damping results in a characteristic 
ratio of convective to conductive heat transfer (the thermal Peclet number, Pe = pc h UR/X, where p is 
uniform density, c h is the specific heat, U is a characteristic velocity and A is the thermal conductivity) is 
small; the dimensionless temperature in this case is then governed by V 2 T=0. Therefore, in this model, 
we assume that convective heat transfer is negligible (i.e., Pe « 1) and that a uniform heat flux density is 
imposed along the vertical wall and constant heat losses through the top and bottom walls of the cylinder 
with the possibility of an additional spatially-periodic heat flux of various amplitude and wavelength 
on the vertical sidewall in some cases. With these assumptions, the energy equation (V 2 T=0) is non- 
dimensionalized by R for length and 2AT for temperature, where A T is the temperature difference between 
the centerline and the vertical (side) wall and the dimensionless analytical solution for temperature is 


T(r,z ) : 


\ 


f r 1 
7 2 

Z 

v 2 , 


+ s 2 I 0 (or)cos(az), 


( 1 ) 


where a = rmHO and n is a known integer, and s 2 are set coefficients, and I 0 is a modified Bessel function 
of the first kind of order zero. The first term on the right hand side of Eq. (1) represents the temperature 
distribution due to a uniform heat flux while the second term represents the temperature distribution due to 
a spatially-periodic heat flux. 


2.2. Flow Problem 

At sufficiently large values of B , the characteristic ratio of the electromagnetic body force to the inertial 
effect (the interaction parameter, N = oB 2 R/(pU), where a is the electrical conductivity of the liquid) is 
very large and, thus, the inertial terms in the Navier-Stokes equations become negligible. In addition 
to the applied magnetic field produced by the magnet, associated electric currents induce an additional 
magnetic field since the materials considered here are electrically conducting. The characteristic ratio of 
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the induced to applied magnetic field strengths is the magnetic Reynolds number, R m = \x pUR, (where 
|l ; is the magnetic permeability of the liquid) and this is very small for crystal growth processes and self- 
diffusion experiments of liquid metals and semiconductors. Thus, neglect of the induced magnetic field 
effect is justifiable in the model presented here. 

Khine and Walker [9] determined the characteristic velocity for magnetically damped axisymmetric 
buoyant convection U - 2pg ft PA77(a.5 2 ), where g (j is the gravitational acceleration (9.81 m/s' 2 ) and [1 is 
the volumetric expansion coefficient for Boussinesq approximation, and this is used to non-dimensionlize 
the velocity. The magnetic flux, reduced pressure (the difference between the total pressure and the 
hydrostatic pressure) with uniform density, the electric current density, and the electric potential are non- 
dimensionalized using B, <jUB 2 R, <jUB, and UBR , respectively. The dimensionless, inertialess equations 
governing the steady, axisymmetric buoyant convection are then 


V*v = 0, (2a) 

V p = Tz+ jx z + Ha' 2 V 2 v, (2b) 

Vy = 0, (2c) 

j = -Vcj) + v xz. (2d) 


Here, Ha = BR(ol\i) m , where p is the viscosity of the liquid, Eq. (2a) is the conservation of mass, Eq. (2b) 
is the Navier-Stokes equation with a buoyancy force and an electromagnetic body force, Eq. (2c) is the 
continuity of electric current density j, and Eq. (2d) is Ohm’s law with an electric potential function (j). We 
eliminated the pressure p, by introducing a stream function, \\ j(r, z ) which satisfies Eq. (2a). 

The boundary conditions are 

5\|/ _ 

V=^"=0, at r = 1 ; 

\|/ = 0, at z = - 10, 10. 

We used the Chebyshev spectral collocation method to solve the equation 

2.3. Diffusion Problem 

The dimensionless axisymmetric governing equation for diffusion is 

5c 

^ — I- Pe m vVc = V 2 c. (4a) 

Here, t is normalized by the diffusion time scale ( R 2 /D ) and c is scaled by the initial concentration C ff Pe m 
= UR/D fj is the characteristic ratio of convective effects to diffusive effects in the process. Here, U is the 
calculated characteristic velocity and differs from case to case. Since there is no mass transfer across any 
boundary, the boundary conditions are 

5c 

qP= 0, at r = 1 ; (4b) 

5c 

■fa = 0, atz = -10, 10. (4c) 

Eq. (4a) is solved for different Pe m with conditions (4b, 4c) using a Chebyshev spectral collocation method 
for spatial discretization as in the flow problem while conventional finite-difference is used for temporal 
discretization. 
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(3a) 

(3b) 

governing \|/(r, z). 


The initial condition for c where tracer forms a thin layer at the bottom of the cylinder is represented by 

c(r,z,t= 0) = exp (-a (z+10) 2 ), (5) 


where a = 3 and dramatically to c ~ 0 around z = -9. The simulated measurements (or output diffusivities 
D) are computed from the resulting time traces of concentration through a straight line fit in the form. 


ln[cj(t)-c 2 (t)] = constant 


( 2 n 
n D 


\ 

t 


( 6 ) 


V 7 

which is known as Codastefano [11] or Harned [12] technique. c 1 andc, are concentration measured atz = 
Z/6 and 5Z/6 along the length of the cylinder, respectively, and the constant depends on the concentration 
profile c(z) at t = 0. 


In the flow problem, 15 collocation points are needed in the radial direction and 40 points in the axial. In 
the diffusion problem, the total number of collocation points in r is 15 and that in z is 80 since the entire 
cylinder is considered. Liquid indium with (3 = 1.02x10 4 K' 1 , p = 6.64x1 0 3 kg-m' 3 and a = 3.02xl0 6 
Sun 1 [13] is used as a model fluid. A self-diffusivity value of D g = 1.48xl0 5 cnF-s' 1 was used as the input. 
Two different Hartmann numbers are considered. For liquid indium these correspond to magnetic field 
strengths of 0.218 T (Ha = 25) and 0.873 T (Ha = 100). The characteristic diffusion time is 1520 s, and 
values of Pe between 0 and 4 are considered here. 

m 

3. Results and Discussion 

The results for five different incoming heat fluxes (see Figs, lb to If) are presented in this section. Case 
I presents self-diffusion with a radial temperature difference only as the driving force in the liquid due to 
the steady, uniform heat flux along the sidewall. Cases II through V present self-diffusion with various 
magnitude of uniform heat flux and superimposed, steady, spatially-periodic heat flux with various 
amplitude and wavelength along the sidewall (i.e., localized hot and cold spots along the sidewall in 
addition to the radial temperature difference in the liquid). The desired incoming heat fluxes for Cases I 
through V are obtained by adjusting s v s 2 and a in Eq. (1). Table I presents summarized results for those 
five different cases. 


Case I: s = 1 and s 2 = 0 in Eq. (1) 

In this case, a uniform heat flux of density 1 that produces a radial temperature difference in the liquid is 
applied through the vertical wall with unifonn heat loss through the top and bottom end walls as in Fig. 
lb. The isotherms show symmetry from z = 0 plane, and are slightly deviated from the horizontal, and 
decrease toward the top and bottom walls. The streamlines for both Ha circulate in counterclockwise 
motion beginning near the vertical wall. 

Figure 2 presents the diffusivity vs. effective Pe m for Ha = 25 and 100. For both Ha, the results lie on one 
curve that bows upward as the effective Pe m increases. From the results, we can see that the convective 
effect is extremely large for a larger effective Pe m . For effective Pe < 2, the predicted convective effect 
in self-diffusivity measurements (i.e., difference in diffusivity from D 0 ) is less than 5% for both Ha = 25 
and 100. The allowable driving force AT, for Ha = 25 is about 0.00373 K while that for Ha = 100 is about 
0.00501 K to prevent convective contamination in excess of 5%. 


45 


Case II: s f = s 2 = 1 and n = 3 in Eq. (1) 

In this case, the uniform heat flux is the same as in Case I which produces a radial temperature difference 
in the liquid while a periodic flux with an amplitude of 0.495 is superimposed at the sidewall to produce 
the localized periodic temperature variation effect. Here, n = 3 represents the number of periodicity in the 
non-unifonn heat flux. The maxima of the periodic flux occur at z = -6.67, 0 (mid-plane), and 6.67 while 
the minima occur at z = -10 (bottom wall), -3.33, 3.33, and 10 (top wall). So, the top and bottom walls, z 
= Z/6 and 5Z/6 are slightly cooler than the rest of the cylinder. 

For Case II, the diffusivity results are plotted in Fig. 3 for both values of Ha. Again, the two curves bow 
upward as the effective Pe m increases. The curve for Ha = 1 00 is slightly higher than that for Ha = 25 . The 
difference is more obvious at a larger effective Pe . For the diffusivity results of less than 5% deviation 
from D n , the allowable AT for Ha = 25 is 0.00373 K and that for Ha = 100 is 0.00548 K. Case II results 
are very close to those of Case I. The spatially-periodic temperature along the sidewall does not affect the 
diffusivities significantly although the presence is observed in the flow pattern. 

Case III: s 1 = 0.01, s 2 = 1, and n = 3 in Eq. (1) 

As shown in Fig. Id, the uniform heat flux with a magnitude of 0.01 (i.e., 100 times smaller than that in 
Case I) and a spatially-periodic flux with an amplitude of 0.495 (same as in Case II) are imposed at the 
sidewall. For Case III, the consequences of the magnetic field for transport can be seen from the diffusivity 
vs. the effective Pe shown in Fig. 3. For both Ha, the curves bow upward as the effective Pe m increases. 
The effect of convection on the measured diffusivity is more prominent at a larger effective Pe m . Also, 
the diffusivity curve for Ha =100 lies above the curve for Ha = 25. For less than 5% deviation from D g , 
the driving AT, for Ha = 25 is about 0.0249 K while for Ha = 100, it is 0. 171 K. Thus, Case III suggests 
that reducing the radial temperature gradient in the liquid by 100 times from its original value (as in Case 
I) has a tremendous effect on the diffusivity results for this model. The effect is more prominent at higher 
Ha for this case (i.e., at stronger magnetic field strength). 

Case IV: s ! = 0.001, s 2 =1, and n = 3 in Eq. (1) 

Here, the unifonn heat flux which produces a radial temperature difference in the liquid is further reduced 
to a magnitude of 0.001 (i.e., 10 times less than that in Case III and 1000 times less than that in Case I) while 
keeping the periodic temperature profile unchanged as in Cases II and III. The diffusivity vs. effective 
Pe m for Case IV is shown in Fig. 3 for both Ha. As in the previous case, the two curves bow upward as 
the effective Pe increases, and the curve for Ha = 100 is higher than that for Ha =25, especially at larger 
effective Pe . The allowable AT in this case for Ha = 25 is 0.0262 K while that for Ha = 100 is 0. 184 K 

m r 

which is a slight improvement from those in Case III. Thus, reducing the magnitude of radial temperature 
gradient by 10 times from Case III does not result in tremendous changes in diffusivity results. 

Case V: s f = 0.01, s 2 =1, and n = 2 in Eq. (1) 

This case is the same as Case III except n = 2 here so that the isotherms decrease and increase two times 
alternatively in between the top and bottom walls (i.e., the periodicity is 2 here). The uniform flux 
which produces the radial temperature difference is 100 time smaller than that in Case I. The maxima of 
spatially-periodic flux occur atz = -10 (bottom wall), 0 (mid-plane), and 10 (top wall) with a magnitude of 
0.207 while the minima occur at z = -5 and 5. So, the top and bottom walls and the mid-plane are slightly 
warmer than the rest of the cylinder. 
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Fig. 1 Model diffusion capillary. Fig. la presents the set up and Figs, lb - If describe the dimensionless 

incoming heat fluxes for Cases I through V. 
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The diffusivities vs. effective Pe m for Case V are plotted in Fig. 3 for both Ha. The curve for Ha =100 is 
higher than that for Ha = 25 and the convective effects are more obvious at larger effective Pe m . In this 
case, the allowable temperature non-uniformity A T for Ha = 25 is 0.0426 K and that for Ha = 100 is 0.222 
K. The results are close to those of Case III with an improvement in sensitivity of allowable temperature 
for the same magnetic field strength. 

4, Conclusions 

The use of an axial magnetic field to offset the detrimental effects of convection on diffusivity 
measurements in liquid metals has been examined using a numerical model. Our goal was to determine 
what magnetic field strengths would be necessary given various temperature non-uniformities to ensure 
that convective transport was less than 5% of the diffusive. Two different flow patterns were observed 
for the various imposed temperature non-uniformities. For a uniform heat flux along the sidewall that 
produces a radial temperature difference in the liquid, the streamlines form a counterclockwise motion. 
For a spatially-periodic heat flux which produces the localized hot and cold spots on the vertical wall, 
the flow is arranged as vertically stacked cells and the number of cells depends on the periodicity of the 
temperature profile due to periodic flux. 
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Fig. 3 Diffusivity vs. effective Pe for Cases II, III, IV and V. 


Table 1 . Summary of results for Cases I - V. 


Case I: uniform heat flux only; d =1, e =0 in Eq. (1) 


\^ = 0.2l8C 
Ha = 25 

<l cS 

V 

1 rmax 1 

(cm/s) 

1 z max 1 

(cm/s) 

U 

(cm/s) 

effective 

Pe 

m 

Dx 1(F 
(cm 2 /s) 


0.00107 

7.34xl0' 6 

6.42x1 0- 5 

9.87xl0- 6 

0.651 

1.49 


0.00373 

2.57xl0 5 

2.25xl0 4 

3.46x10 s 

2.28 

1.55 


0.0107 

7.34xl0' 5 

6.42x10 4 

9.87xl0' 5 

6.51 

1.99 


0.0426 

2.94x1 0- 4 

2.57xl0- 3 

3.95xl0- 4 

26.04 

9.57 

= 0.8737' 
Ha = 100 

AT 

(K) 

\ v 

1 rmax 1 

(cm/s) 

1 z max 1 

(cm/s) 

U (cm/s) 

effective 

Pe 

m 

DxlO 3 

(cm 2 /s) 


0.00171 

2.85xl0' 6 

7.48xl0- 5 

9.87xl0' 7 

0.758 

1.49 


0.00501 

8.35xl0 6 

2.19xl0 4 

2.89xl0 6 

2.22 

1.55 


0.0171 

2.85xl0' 5 

7.48xl0- 4 

9.87xl0- 6 

7.58 

2.11 


0.0854 

1.43x1 O' 4 

3. 74x1 0' 3 

4.93xl0' 5 

37.9 

16.0 
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Case II: uniform flux ~ periodic flux; d= 0.01, e = 1 and n = 3 in Eq. (1) 


= 0.2i8C 

AT 

]v rnv 

V Tm 

U 

elfective 

Dx lO 3 

Ha = 25 

(4 

(cm/s) 

(cm/s) 

(cm/s) 

Pe 

m 

(cm 2 /s) 


0.00533 

8.78xl0' 6 

6.41xl0- 5 

4.93xl0' 5 

0.652 

1.49 


0.0249 

4.11x10 s 

3.00x10 4 

2.31xl0 4 

3.04 

1.55 


0.0639 

1.05xl0' 4 

7.70x1 0- 4 

5. 92x1 O' 4 

7.82 

1.84 

^^5 = 0.8737 

AT, 

\ v 

\ v 

U 

effective 

£>xl0 5 

Ha = 100 

(K) 

(cm/s) 

(cm/s) 

(cm/s) 

Pe 

m 

(cm 2 /s) 


0.0427 

8.46x1 O' 6 

6.04xl0' 5 

2.47xl0' 5 

0.612 

1.49 


0.171 

3.38x10 s 

2.42x10 4 

9.87x10 s 

2.45 

1.55 


0.684 

1.35xl0- 4 

9.67x1 0- 4 

3.95xl0' 4 

9.79 

2.2 

Case III: uniform flux ~ 100 * periodic flux; d = 

e = 1 and n = 

= 3 in Eq. (1) 

= 0.218 T 

AT 

V 


U 

effective 

Z)xl0 3 

Ha = 73 

(JO 

(cm/s) 

(cm/s) 

(cm/s) 

Pe 

m 

(cm 2 /s) 


0.00107 

5. 92x1 O' 6 

7.56x1 0- 5 

9.87xl0- 6 

0.766 

1.49 


0.00373 

2.07x10 s 

2.65xl0 4 

3.45x10 s 

2.68 

1.55 


0.0107 

5. 92x1 0' 5 

7.56xl0- 4 

9.87xl0' 5 

7.66 

1.99 

= 0.8737' 

AT 

|v r 

\v 

U 

elfective 

27x1 0 3 

Ha= 100 

(JO 

(cm/s) 

(cm/s) 

(cm/s) 

Pe 

m 

(cm 2 /s) 


0.00171 

2.83xl0' 6 

7.64xl0- 5 

9.87xl0' 7 

0.774 

1.49 


0.00548 

9.07xl0 6 

2.45x10 4 

3.16xl0 6 

2.48 

1.55 


0.0171 

2.83xl0' 5 

7.64x1 0- 4 

9.87xl0' 6 

7.74 

2.12 

Case IV: 10 * uniform flux ~ periodic flux; d = 0.001, e = 1 and 

n = 3 in Ec 

•(1) 

\77 = 0.2187’ 

AT 

\ v 

\ v 

U 

effective 

27x1 0 3 

Ha = 33 

(JO 

(cm/s) 

(cm/s) 

(cm/s) 

Pe 

m 

(cm 2 /s) 


0.00533 

8.78xl0' 6 

6.17xl0- 5 

4.93xl0' 5 

0.624 

1.49 


0.0262 

4.31x10 s 

3.03x10 4 

2.42xl0 4 

3.07 

1.55 


0.0639 

1.05xl0- 4 

7.40x1 0- 4 

5. 92x1 O' 4 

7.50 

1.84 

^^77 = 0.8737’ 

AT 

|v m . 


U 

effective 

27x1 0 3 

Ha = 100 

(JO 

(cm/s) 

(cm/s) 

(cm/s) 

Pe 

m 

(cm 2 /s) 


0.0427 

8.61xl0' 6 

5. 06x1 O' 5 

2.47xl0' 5 

0.513 

1.49 


0.184 

3.71x10 s 

2.18xl0 4 

1.06x10"* 

2.21 

1.55 


0.513 

1.03xl0' 4 

6.07xl0- 4 

2.96x1 O’ 4 

6.15 

1.87 


Case V: uniform flux ~ periodic flux; d= 0.01, e = 1 and n = 2 in Eq. (1) 


\^ = 0.2l8C 
Ha = 75 

AT 

r 

(K) 

V 

1 rmax 1 

(cm/s) 

v 

1 z max 1 

(cm/s) 

U 

(cm/s) 

effective 

Pe 

m 

Dx 1(P 
(cm 2 ls) 


0.0107 

1.04x1 0' 5 

8.39x10 s 

9.87xl0' 5 

0.850 

1.49 


0.0426 

4.14x10 s 

3.35xl0 4 

3.95xl0 4 

3.40 

1.55 


0.107 

1.04x1 O' 4 

8.39xl0- 4 

9.87xl0- 4 

8.50 

1.86 

= 0.8737' 

AT 

r 

~~ [v 

1 rmax 1 

|v 

1 zmax 1 

U 

effective 

27x1 0 3 

S 

II 

© / 
© / 

(K) 

(cm/s) 

(cm/s) 

(cm/s) 

Pe 

m 

(cm 2 /s) 


0.0854 

l.OOxlO' 5 

6.91xl0- 5 

4.93xl0- 5 

0.700 

1.49 


0.222 

2.60x10 s 

1.80x10 4 

1.28xl0 4 

1.82 

1.55 


0.854 

1.00x1 O' 4 

6.91xl0- 4 

4.93xl0- 4 

7.00 

2.26 


In all cases, the value of the simulated diffusivity increased with increasing effective Pe (i.e., with 
increasing temperature difference, and, thus, convective velocity magnitude). The dependence of the 
transport conditions on the nature of the thermal boundary conditions was such that for a unifonn heat 
flux along the sidewall which produces a radial temperature difference in the liquid, the diffusivity vs. Pe m 
were indistinguishable for //a = 25 and 100. The diflfusivities for Ha =100 were higher than those for Ha 
= 25 in the presence of a periodic temperature profile produced by a periodic heat flux superimposed on a 
uniform heat flux (Cases II, III, IV, and V). 

Case I which is uniform heat flux only (i.e., the radial temperature difference in the liquid is the only 
driving force) seems to be the most sensitive while Case V is the least sensitive among five different 
conditions considered in this model. From the results, one can conclude that a particular temperature 
profile at the sidewall (produced by a combination of unifonn and non-unifonn heat fluxes) is required to 
result in desired diffusivities (i.e., to operate within the allowable temperature non-uniformities) in self- 
diffusion experiments for this numerical model. For a fixed value of effective Pe , A T increases with 
increasing Ha (the magnetic field strength). Thus, stronger magnetic fields can tolerate a stronger driving 
force A T in self-diffusivity measurements for same convective conditions. 
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